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Abstract 

The production of hard photons and neutral pions in 190 MeV proton induced 
reactions on C, Ca, Ni, and W targets has been for the first time concurrently 
studied. Angular distributions and energy spectra up to the kinematical limit are 
discussed and the production cross-sections are presented. From the target mass 
dependence of the cross-sections the propagation of pions through nuclear matter is 
analyzed and the production mechanisms of hard photons and primordial pions are 
derived. It is found that the production of subthreshold particles proceeds mainly 
through first chance nucleon-nucleon collisions. For the most energetic particles the 
mass scaling evidences the effect of multiple collisions. 
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1 Introduction 



Forming hot and compressed nuclear matter by colliding heavy ions provides 
the today only known experimental mean to explore the phase diagram of nu- 
clear matter. During the collision the longitudinal momentum of the projectile 
is converted by multiple nucleon-nucleon collisions into transverse momentum 
as well as through the production of secondary particles, leading to the cre- 
ation of a hot and compressed reaction zone. At projectile energies below 
the free pion production threshold (£^° } =280 MeV) pions and hard pho- 
tons {E 1 > 30 MeV) carry information on the initial compressional and high 
density phase of the heavy-ion collision which might evolve towards a hot 
thermalized nuclear system [1,2]. 

In first order, subthreshold pions and photons are created following either sin- 
gle or several nucleon-nucleon collisions, the nuclear medium providing the 
additional energy needed to surpass the production threshold. Therefore sub- 
threshold pions and photons can be viewed as a snapshot of the dynamical 
phase-space occupancy of participant nucleons at the early instant of the col- 
lision. However, extracting the information relevant to the thermodynamical 
state of the system is far from being straightforward. One must rely on mod- 
els which describe the complicated many-body problem of a heavy-ion col- 
lision. Alternatively, a purely empirical approach consists in the comparison 
of various observables between heavy-ion collisions and proton-nucleus colli- 
sions. Since a compressional phase is very unlikely to occur in the latter, the 
in-medium nucleon-nucleon collisions, i.e. the longitudinal momentum dissi- 
pation, can be studied without the complication due to the dynamics and 
thermodynamics of heavy-ion collisions. 

For the present measurement the proton beam energy was chosen to be 190 
MeV, resulting from a compromise between the concurrent requirements to 
accumulate a usable set of data and to be well below the pion production 
threshold. The energy spectra, up to the maximum energy available in the 
reaction for the production of a single particle, were measured together with 
their angular distributions. The target mass dependence of the production 
cross sections, studied as a function of the energy of the produced particle, 
is exploited to characterize the production mechanism and, in the case of pi- 
ons, their propagation properties through nuclear matter. It is concluded that 
the production mechanism of subthreshold particles changes with increasing 
particle energy from a production in individual nucleon-nucleon collisions to 
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mechanisms involving more than two nucleons. Although several similar stud- 
ies have been published earlier [3-8], the information collected in the present 
work is unique as it combines for the first time the measurement of com- 
plete energy spectra with the concurrent observation of strongly interacting 
particles, pions, and non-interacting particles, photons. 



2 Experimental setup 

The 190 MeV proton beam was delivered by the superconducting cyclotron 
AGOR at KVI. Natural C, Ca, Ni and W targets, respectively 18.1, 11.2, 6.7 
and 5.0 mg/cm 2 thick, were irradiated with 10 to 25 nA beams in bunches of 
2 ns at a rate of 60 MHz. Within these conditions the interaction rate was 
respectively for the four targets 0.23, 0.13, 0.09, and 0.05 nuclear interactions 
per beam pulse. These values set the counting rate in an individual photon 
detection module to 1 kHz and induced negligible double hits and random 
coincidences. Direct photons and decay photons from neutral pions (ir° — > 77, 
BR=98.8%) were detected with the photon spectrometer TAPS configured in 
six blocks of 64 BaF 2 modules, each module being associated with a charged- 
particle-veto detector. The blocks were positioned at the most forward angles, 
allowed by technical constraints, on both sides of the beam. They covered 
a quasi-continuous polar angle range above 9 = 60° and an azimuthal angle 
range of ±20°, hence covering approximately 15% of 4-7T. The distance from the 
target to the active front face of the blocks (d = 66 cm) offered a flight path 
large enough to discriminate photon hits against baryon (mainly neutron and 
proton) hits through their time-of-flight measurement. The energy deposited in 
the BaF 2 modules was calibrated and the gain continuously monitored during 
the experiment using the 38.5 MeV peak energy-deposition of cosmic muons. 
Time-of-flight was calibrated using the prompt photon signal and the known 
time interval between proton beam pulses (16.7 ns). To enrich the recorded 
events with photons and neutral pions the trigger electronics was set up to 
select events in which at least one, respectively two, neutral hits (defined as a 
hit in a BaF 2 module and no hit in the charged-particle-veto detector) were 
present, depositing at least 15 MeV equivalent photon energy. A downscaled 
trigger consisting of a neutral hit of at least 0.4 MeV served as minimum- 
bias trigger. An exhaustive description of the TAPS detector, its associated 
electronics and trigger scheme can be found in [2]. 

In the offline analysis photon hits were identified by combining the informa- 
tion delivered by the BaF 2 modules and by exploiting the characteristic pulse 
shape and time of flight of photons. The direction and energy of the im- 
pinging photon was derived from the analysis of the electromagnetic-shower 
topology[3,4] which provided also additional discrimination criteria against 
baryons and cosmic muons. Neutral pions were identified, with a full width at 
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half maximum of 10%, through an invariant mass analysis of identified photon 
pairs. In the resulting invariant mass spectrum combinatorial background can 
be safely neglected because of the weak 7 multiplicity at the considered bom- 
barding energies (M 7 ss 1CT 3 ). Photon pairs with an invariant mass between 
100 and 160 MeV were adopted as stemming from the decay of the produced 
neutral pions. Close to one million pions were detected and identified for each 
target. 

For every identified neutral pion, the 4-momentum was calculated through 
a constrained minimization method using the 4-momentum of the two pho- 
tons. This special energy reconstruction method was developed[5] because the 
standard technique [6] leads to unphysical results [7] (pion energies above the 
kinematical limit) in our case, where the pions i) are of very low energy and 
ii) take away a large (up to 100%) part of the total energy available in the 
proton+nucleus system. Performances of this new method have been carefully 
checked using GEANT simulations [8] of mono-energetic neutral pions. Unlike 
the standard method, our minimization method leads to a non-symmetric re- 
sponse function (Fig. 1), with a reconstructed mean energy which is about 
88% of the true value, and an energy resolution depending on pion energy 



The neutral pion efficiency was calculated with the help of GEANT simula- 
tions taking as input an extrapolation D(K w o, £\o) to the full solid angle of the 
measured double differential cross-sections d 2 a/d£l n odK n o. D^K^o,^}^)^ = 
da/dK^o (angular distribution has been chosen isotropic, in agreement with 
the data) where da/dK w o was obtained through the following iterative proce- 
dure. The distribution of the measured kinetic energy, K n o, was parametrized 
as : 



Expression 1 was selected as the initial guess of the iteration procedure and 
fed into the simulation. The parameters a,T 1 ,T 2 , K were then tuned until 
the output of the simulation reproduces the measured distribution. 

The photon efficiency was calculated using the same scheme, with the following 
input distribution : 



(Fig. 1). 




(1) 




The resulting pion and photon efficiencies are reported in Table 1. 
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3 Cross-sections 



The neutral pion total cross sections were calculated as: 



o. 




C D x N r J w /e. 



where N™ w is the raw number of identified pions, Cp is a normalization factor 
including target thickness, beam intensities and trigger conditions, and e 4 ! is 
the global efficiency. The resulting cross-sections <r 4 o are reported in Table 1. 

To obtain direct photon cross-sections an additional treatment is required, as 
the contribution of the decay photons from neutral pions must be subtracted 
from the total photon cross-section. The latter is obtained as 



where e 47r = 12% is the efficiency for photons and N™ w is the raw num- 
ber of identified photons. The contribution of the decay photons was esti- 
mated with the help of GEANT simulations taking as input the extrapolated 
double-differential cross-section D(K n o : Q n o). The contribution of decay pho- 
tons amounts to about 30% of the photon spectrum (Fig. 3). 

From the photon cross section the probability P 7 to produce a single pho- 
ton in an individual proton-neutron collision was evaluated according to the 
prescription of [9] . We observe a good agreement of this probability with the 
systematics presented in [9] (averaged for Ca, Ni, W targets, the P 7 equals to 
(7.4 ± 1.4) x 10~ 4 , compared to 9.4 x 10~ 4 expected from the systematics) with 
the exception of the C target (factor 2 below the systematics). The neutral 
pion cross sections can be compared to the results of [10] obtained with proton 
beams of 200 MeV kinetic energy. While we observe the same (within error 
bars) mass scaling law of the total cross section for neutral pion production, 
the cross section values reported in [10] are higher by a factor of approxi- 
mately 3. Such a strong rise can not be attributed to the 10 MeV difference 
in beam energies. The problem of inconsistencies in the cross section for pion 
production in proton-nucleus reactions is discussed in detail in [11]. 

Thanks to the good statistics accumulated, errors in this experiment are 
mainly systematic in nature, and amount to 10% for Cd, 7% for e 4 o, 20% 
for e 471 ", leading to total relative errors of about 32% for the direct photon 
cross sections and 12% for the total pion cross sections. 

The energy spectra of particles (Figs. 2 and 3) produced from the four different 
targets are conveniently compared if one defines a reduced energy, e r , which 
represents the fraction of the available kinetic-energy, K max , carried away by 
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the particle. This energy is calculated (Table 2) assuming that the reaction 
proceeds through the fusion of the proton and the target nucleus: K max = 
y/s — (A + 1) M N — Q, where \/s represents the total energy available in the 
center-of-mass, M N the nucleon mass and Q the fusion reaction Q-value. For 
photons K = E 1 . Within this representation (Figs. 2 and 3, right panels) 
particles with e r — 1 take away all the available energy and the mechanism 
would be pionic fusion [12] for pions or radiative capture for photons. 



4 Analysis of target mass dependence 

The dependence with the target mass A of the pion and photon cross-sections 
can be described by a power law, da/dx oc A a ^ , where x is either the po- 
lar angle of the emitted particle or its reduced kinetic energy. This widely 
used parametrization at relativistic [13-16] as well as ultrarelativistic [17,18] 
energies for a variety of produced particles gives valuable indications on the 
particle-production mechanism. Since one is interested in the primordial par- 
ticle production, one must first understand how the primordial production 
is modified by subsequent absorption or re-scattering during the propagation 
through the nuclear medium. This has been done by examining the target 
mass dependence of the differential cross-section, da/d6 n o, in terms of an el- 
ementary geometrical model which tracks the pion through the nucleus after 
its creation [5]. Within this picture our data are best described if one as- 
sumes that: i) primordial pions are produced in a surface shell of the nucleus 
with thickness \ p = 2.4 fm, which is equal to the proton mean free path, 
and ii) the pion mean free path, independent of the pion energy, is equal to 
\ w = 6 fm in good agreement with the values established elsewhere [19]. The 
first assumption implies that pions detected in the backward hemisphere are 
quasi unaffected by re-interaction in the medium and those pions can thus be 
considered as primordial pions. 

The target mass A dependence of neutral pion and photon cross-sections was 
analyzed in terms of the power-law parameter a. The carbon data have been 
excluded from this analysis because of its extraordinary behavior (see also 
[20]). All neutral pions and photons with energy above 92 MeV were pro- 
duced below the free NN threshold. For hard photons produced below thresh- 
old the a parameter (Fig. 4) continously increases from a = 2/3 to values 
close to or even above one. This behavior is interpreted as an evolution of 
the particle-production mechanism from a surface production for the least 
energetic particles towards a volume production for the most energetic parti- 
cles. In other words, the observed a dependence witnesses a transition from 
particles produced in first-chance nucleon-nucleon collisions to particles pro- 
duced in multiple successive binary collisions. A possible mechanism for the 
latter is the process of photoabsorption of subthreshold pions in the nuclear 
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medium (NN — > NNtt and subsequently Ntt — > ^7), as proposed in [21]. 
The a dependence with the photon energy (Fig.4) indicates another strong 
rise when the energy of photons decreases. This can be understood invoking 
again the participation of all nucleons with the difference that in this case 
there is enough energy available, not only in first chance collisions but also in 
secondary collisions. A detailed analysis of this photon energy region was per- 
formed by another TAPS experiment during the same experimental campaign 
at KVI [20]. 

In the case of neutral pions such firm conclusions can not be drawn from the 
available experimental data. Within error bars, a rise of the a parameter to- 
wards unity or a constant value of a = 2/3 can not be distinguished. It must 
be emphasized that, compared to photons, pions are less adequate probes for 
a model-independent analysis. Their production often involves baryonic reso- 
nances as intermediate steps [22], and the energy dependence of the power-law 
parameter a might be affected by the energy variation of the pion absorption 
length, an effect not accounted for in the simple geometrical model used here. 



5 Conclusion 

A complete data set on the production of neutral pions and photons in re- 
actions of 190 MeV protons with targets ranging in mass from A = 12 to 
A = 184 has been obtained. Energy spectra up to the kinematical limit have 
been presented. A simple geometrical model, assuming neutral pion produc- 
tion in a surface shell of the target nucleus and subsequent pion propagation 
with the mean free path equal to = 6 fm, agrees well with the observed 
angular distribution of neutral pions. The analysis of the mass scaling law of 
subthreshold particle cross-sections supports the mechanism of particle pro- 
duction in first-chance nucleon-nucleon collisions, except in the case of most 
energetic particles (up to the kinematical limit), where the role of multi-step 
processes becomes evident. 

The present data furnish new and valuable information for tuning unknown pa- 
rameters of dynamical models for heavy-ion collisions[l], like in-medium cross 
sections or the implementation of the Pauli principle. Further evidence for new 
production mechanisms like pn — > dn (7) or nN — > N'j [21] has been given. 
The complicated and hardly understood production mechanism of hadronic 
particles, like pions and kaons [23], affected by reabsorption and rescatter- 
ing, stresses the importance of hard-photon measurements as presented in 
this letter. Once hadronic and electromagnetic probes are well understood, 
the dynamical transport model codes may provide a better description of the 
particle production in heavy-ion collisions, which is the base to unravel the 
thermodynamic properties of nuclear matter. 
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Table 1 

Cross sections of neutral pions and photons measured in p+A reactions at 190 
MeV. N™ w is the number of neutral pions recorded, Co takes into account target 
thickness, beam intensity and trigger rates, cr^o = Co x N™ w /e%l is the total neutral 
pion cross-section, calculated using estimated neutral pion global efficiency. The 
direct photon cross-section a^ irect = a l ° tal — £j^°^ 77 is then deduced from the total 
photon cross-section by subtracting contribution from decay photons. All photon 
cross-sections are given for E l °- b > 40 MeV. 





C 


Ca 


Ni 


W 


C D x 10 6 (/ib) 


1.45 ±0.15 


10.2 ± 1.0 


2.38 ±0.24 


12.5 ±1.2 


N™ w x 10" 5 


3.042 ± 0.005 


2.017 ±0.004 


9.870 ± 0.01 


4.166 ± 0.006 


45 (%) 


2.65 ±0.18 


2.43 ±0.18 


2.43 ±0.18 


2.44 ±0.18 


af (fib) 


16 ±2 


85 ± 10 


97 ± 12 


215 ±26 


(%) 


0.12 ±0.02 


0.12 ±0.02 


0.12 ±0.02 


0.12 ±0.02 


a l ° tal (/xb) 


88 ±20 


520 ± 120 


567 ± 130 


1383 ± 320 


(fib) 


28 ±3 


149 ± 18 


170 ± 20 


376 ± 45 


a direct ^ 


60 ±20 


372 ± 120 


397 ± 130 


1007 ±310 



Table 2 

Reaction Q-values (from [24]), used to compute the maximum available K max in 
A(p,X) reactions : Q gg = AM A+1 - AM A - AMi + m x where X = tt°, 7. 



Reaction 


Q gg (MeV) 


K max (MeV) 


i 2 C 3 . s .(p,vr )fN 9 . s . 


133.0 


41 


2oCa g . s . (p, vr )2{Sc 9 . s . 


133.9 


51 


lNi g . s .(p,7T )|Cu g . s . 


131.5 


55 


ifW s . s .(p,vr°)i 5 Re 9 .3. 


129.6 


59 


i 2 C,. s (p, 7 )PN,. s . 


-1.94 


176 


!gCa s . s .(p, 7 )fiSe 5 ., 


-1.08 


186 


|Ni,. s .(p,7)|C%. s . 


-3.42 


190 


^ 4 W,. s .(p,7)^l 5 Re 9 . s . 


-5.40 


194 
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Fig. 1. Left: typical simulated response function of TAPS to mono-energetic neutral 
pions of 25 MeV, represented by the difference between the true initial pion kinetic 
energy and the reconstructed pion kinetic energy modified by the detection system 
and the reconstruction altogether. Right: Performance of the neutral pion momen- 
tum reconstruction method, as a function of the neutral pion kinetic energy (top: 
mean of the reconstructed energy, bottom: resolution). 
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Fig. 2. Neutral pion energy spectra measured in p+A reactions at 190 MeV, as 
a function of kinetic pion energy (left) and of the reduced neutral pion energy 
e w = K^/K max (right, see text). 
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Fig. 3. Left: contribution of decay photons from ir° (dotted line) in the total (thin 
solid line) photon spectra, and the resulting direct photon spectra (bold solid line). 
Right: direct photon energy spectra measured in p+A reactions at 190 MeV, as a 
function of the reduced photon energy (from bottom to top: C, Ca, Ni and W). 
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Fig. 4. Evolution of the a parameter as a function of the particle energy, for pho- 
tons (open symbols) and neutral pions (closed symbols). The arrow indicates the 
separation between below- and above-threshold photons. 
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